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Abstract 

An unified systematic of eiastic (anti)proton-proton scattering data is 
proposed based on a simple expression for the process amplitude - 
f(q) = A [{ndq) / sh{TTdq)][i Ji{Rq)/{Rq) + pJo{Rq) + ...]. The parame- 
ters R and d are obtained at a fitting of (anti)proton-proton experimental data 
on differential cross sections from Piab ~ 1 GeV/c up to Tevatron energies. 
The fitting gives extra-ordinary good results, /NoF ~ 1 of below at \t\ < 
^SJ ■ 1.75 (GeV/c)^. An extrapolation of the parameter's energy dependencies to 

the LHC energies allows a good description of the Totem data up to the second 
. diffraction maximum. Predictions for other LHC energies are presented also. 

The amplitude provides one with parameterizations of total and elastic 
cross sections. Its impact parameter representation corresponds to the 
2-dimensional Fermi- function - 1/[1 -I- exp{{b — R)/d)], which is very useful 
^\ ' for Glauber calculations of nucleus-nucleus cross sections at super high energies. 

(N" 

It is shown for the first time that experimental high \t\ elastic scatter- 
(— I . ing data have a weak energy dependency. This allows to describe high \t\ tail 

^ I of the Totem data. 

p : 

(D . Introduction 

^; 

According to the experimental data on elastic (anti)proton-proton scattering data in the Coulomb-nuclear 
ly-^ . interference region, the nuclear elastic scattering amplitude in the momentum representation, f{q), has 

^ ■ a small real part and a large imaginary one. Correspondingly, the amplitude in the impact parameter 

■"sj" . representation, ^{b), has a large real part and a small imaginary one. The amplitudes are connected by 

9P " the Fourier-Bessel transform: 

ON; 

^ ■ ^('^^ = i y '^i{b)dH = I b MS) i{b)db, 

7(6) e^^f{q)d^q = - T q MS) fiq)dq, (2) 

47r /m/(0) = 47r / b Re j{b)db, (3) 
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cr"' = 2n b |7(6)pd5, (4) 
Jo 

d'^a/dq^ = \f{q)\^, Ref{0)^0 Imj{b) ^ 0, (5) 

where q is the momentum transfer, the corresponding 4-momentum transfer - t = —(f'. 

Any function existing on the semi- infinite interval, [0 — oo] , can be represented in a first approximation 

as 

7(6) ~ai [e(i?i-6)-e(0-6)], 

where 



Ri^2 6^ j{b)db/ / b 7(6)d6, 
Jq Jo 



oo 
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ai = (2/i?2) / b -i{b)db. 
Jo 

In a next approximation it is needed to consider 2 functions: 

7(Li)(&)=7W-0(i?i-6), &<^i, 

l(Ri){b) ^lib), b>Ri. 
The function 7(fli)(6) can be approximated as, 

7(iii)W = a2 [e(i?2-6)-e(i?i-6)], 



,2 ..T /r r,2 -fl 



R2 = 2I3/I1 - R^, 02 

The procedure can be repeated one more giving new functions - 7(L2)(^) a-nd 7(^2) (^)- One can obtain 
the results continuing the procedure, 

00 00 
7(6) = ^ a, [e(i?, - 6) - e(i?,_i -b)] +J2 7(L.) ^0 = 0. (6) 

i=l i=l 

The corresponding f{q) is (see [1]), 

Each term in the second series in Eq. [5] can be decomposed on the wavelets |1] , especially, on the 
Haar wavelets. This will lead to an appearance of new 0-functions. Each term of the series in Eq. [7] can 
be decomposed in the Talor series on {Ri — in the vicinity of Ri. Because derivatives of the Bessel 
functions, Jq and Ji, are expressed through each other, the final expression will be: 

/('Z) = ,/i(?)^^^^ + f2{q)MRiq) + (8) 

where /i, /2, are non-oscillating functions. 

The first term of the Eq. [8] was considered [3l|4] many years ago in the Strong Absorption Model (SAM) 
in application to hadron-nucleus and nucleus- nucleus scattering at low energies. Various form of the 
smearing function, /i, were proposed in that time without solid physical foundation. fi{q) = ndq/ sh{TTdq) 
was used in papers [1] by W.E. Frahn et al. Recently, this form of the smearing function was obtained 
[S] at the Fourier-Bessel transform of the symmetrized 2-dimensional Fermi-function, 

7(6) ^ ^ 



fsFpiq) sinh{ndq) Rq ^ 2q'^ sinh(ndq) ^tanh(ndq) ^ o( ^" ( ) 

The symmetrized function is very close to the ordinary Fermi-function, 

where we introduce the coefficient, A, to have more general expression. Usually it is assumed that A ^ \. 
The function (fTU)) was used in the paper [5] by P. Brogueira and J. Dias de Deus for a description of elastic 
pp-data at ^/s — 14, 20, 53 GeV, and pp-data at y/s = 546, 630, 1800 GeV. "Unexpected qualitative 
agreement with the data was found" by the authors. 
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Because the expression [8] is general, we accept as a working hypothesis that the elastic scattering 
amplitude can be described by the expression, 



Rq 2 \tanh{ndq) ) [Rqf 



sinh{'Kdq) 

+ TT^dP/i) {p - ai{-Kdqf) Jo{Rq)/2} 



(11) 



The parameter A is a normalization constant. We consider it as a free parameter in order to take into 
account possible normalization error of experimental data. 

/9 is a ratio of the real to imaginary parts of the elastic amplitude at zero momentum transfer. The 
functional form of the real part, Jo{Rq), is obtained with a help of the derivative dispersion relations (see 
[7]) applying them to the Eq. IHlat 9 — 0. The relations give also other terms which go to as {ndq)'^ at 
q 0. Thus we include them efficiently having the last term in Eq. [11] 

The parameter R is responsible for a position of the first difFractional minimum. The parameter d 
determines the slope of the differential cross section. A filling of the diffractional dip is connected with 
the parameters p and ai. p parameter is not enough to do this. Thus, we introduce a dependence of the 
ratio of the real to imaginary parts of the amplitude on the momentum transfer - q. 

Below we present in Sec. 1 results of a fitting which we make to find the parameters R, d and oi. We 
have used here a lot of experimental data on the differential cross sections of elastic pp- and pp-scatterings. 
It is found that the parameters have rather simple energy dependencies. 

The dependencies are applied in Sec. 2 at an extension to the LHC energies. The extrapolated values 
of R, d and ai together with the expression [TT] allow to describe the Totem data. According to our 
estimation (see Sec. 1) the parameterization is valid at \t\ < 1.75 (GeV/c)^, thus the data are described 
up to the second diffraction maximum. 

In Sec. 3 we consider a behaviour of the differential cross sections of elastic scattering at large \t\. 
Here we present a collection of various experimental data which shows that the cross sections have a weak 
energy dependencies. We propose a simple parameterization of the high |<| tail of the cross sections. This 
allows to describe the Totem data at large \t\. 

At last, in Sec. 4 we consider possible applications of the proposed approach. As seems to us, the 
most important one is its application in calculation of the Glauber cross sections of hadron-nucleus and 
nucleus- nucleus interactions needed for cosmic ray studies. Usually, it is assumed in the calculations, 
that 7(5) has the gaussian form. It would be well to use the Eq. (TO] 

An another application is a usage of the Eq. \W\ in a calculation of quark-gluon string multiplicity 
distributions in hadron-nucleon, hadron-nucleus and nucleus-nucleus interactions within the framework 
of the Quark-Gluon-String Model (QGSM) JW . The analogous application can be foreseen in the high 
energy model like HIJING 9J combined soft and hard interactions. 

The list of the possible applications is not completed. We hope it will be extended in the future. 



1 Fitting of experimental data 

It is obvious that the proposed expression 1111 cannot be applied at all q because cross sections predicted 
by it decrease exponentially with a q growth. At the same time it is known that experimental data show 
a fall down as l/|i|^~* where, as it is expected, hard QCD-processes dominate (|i| — q^). Thus we have 
to restrict the application region of our approach. To find a "border" between hard and soft interactions 
we have undertook a research, results of which are presented in Fig. 1. 

The solid (black) line there shows a fitting results without restriction on \t\ with 5 free parameters 
[A, R, d, p and ai). As seen, the fitting reproduces the data at large |t|, but underestimates the data 
rather strongly at small \t\. At the maximum allowed \t\ = 2 (GeV/c)^ we have the result presented by 
the dashed (green) line. In the case, the fitting curve starts to deviate from the data only at \t\ = 1.75 
(GeV/c)^. At \t\max = 1-5 (GeV/c)^ there is no problems (see dotted (red) line). Thus, as a compromise, 
we have estimated \t\max as 1.75 (GeV/c)^. 

In the following we use a lot of experimental data from 1 GeV up to the Tevatron energies. Some 
part of them were taken from data-base [TT] created by J.R. Cudell, A. Lengyel and E. Martynov. It was 
described in the paper |12| . A complete list of references is given in the Appendix. Because a correlation 
between the parameters was very strong at a fitting of the small scattering angle region, we selected data 
where the dip region was presented. This restricted the set of the experimental data. 
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O 1 2 3 

|t| (GeV^) 

Figure 1: Differential cross section of the elastic pp- scattering at -y/s = 62.3 GeV. The points are the 
experimental data ^Oj. Lines are the fitting results. 



A good results of the fitting were obtained for pp- interactions [TB] ■ More complicated situation took 
place with the fitting of the pp-data. In order to reduce the number of the free parameters, we have 
fixed p using the following parameterizations of the corresponding experimental data from the PDG 
data-base [Tl] : 

o A on 

Ppp = 0.135-^-(---t-2^. (12) 



Ppp — 0.135 — 



2.26 



(13) 



We assume that this allows us to attract indirectly an additional experimental information because p 
values were measured in the independent experiments - in the Coulomb-nuclear interference region at very 
small q. In principle, p can be obtained at the fitting using very exact measurements in the dip region. 
Such data were presented by the EDDA collaboration [TS] for projectile proton momenta, Piab^ from 1.1 
GeV/c up to 3.3 GeV/c. A clear change of the slope of the experimental cross sections was observed at 
Piab > 2 GeV/c and 9cms > 60°. Our obtained p values for the data at Piab > 2 GeV/c are in a reasonable 
agreement with other experimental data [Tl]. Below 2 GeV/c the 4-parameter fit was unstable due to 
the parameter's correlations. Thus we estimate the lower energy boundary of the application region of 
the present approach as 2 GeV/c for pp-interactions. For pp-interactions the boundary can be smaller. 
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Figure 2: Quality of the fitting of pp and pp experimental data. 



4 



With all the above given restrictions we have extra-ordinary good fit. /NoF = 1156/1489 ~ 0.78 
for pp-interactions at Piab > 9.9 GeV/c, and x^/NoF = 856/675 ~ 1.27 for pp-ones at Piab > 8 GeV/c. 
Thus, it seems to us, that we can say about the unified systematic of all high energy (anti)baryon-baryon 
elastic scattering data. 

A quality of the fit is shown in Fig. 2. As seen, most of the pp-interactions data are described quite 
well. The situation is more complicated for the pp-data especially at low energies. At high energies, the 
quality of the fitting of pp-interaction data becomes better. 

Some fitting results in a comparison with experimental data are presented in Fig. 3, 4. We show there 
the experimental data at all measured values of |i| and our results extended outside the fitting region 
{\t\ < 1.75 (GeV/c)^) in order to demonstrate a necessary to include a description of the large angle 
scattering. As seen, we reproduce the cross sections for pp-interactions in the fitted region of \t\. The dip 
position is reproduced also. 




Figure 3: Differential cross sections of elastic proton-proton scatterings. The points are the experimental 
data (a - [TB], b - [T7], c - [TU], d - [H]). The solid lines are the fitting results without the restriction on 
|t|. 




Figure 4: Differential cross sections of elastic antiproton-proton scatterings. The points are the experi- 
mental data??? (a - [Tl],b - [20], c - [H], d - [H]). The solid fines are the fitting results without the 
restriction on \t\. 

One can see that the diffraction minimum in pp-interactions connected with the first zero of Ji{x) 
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shifts to low values of |i| with energy growth. This signals that R is an increasing function of s. The 
filling of the dip is caused by a variation of p. At Piab > 1-5 GeV/c and Piab < 200 GeV/c p is negative. 
At Piab ^ 200 GeV/c p is closed to zero. At larger energy it is positive. So, the fiUing of the dip depend 
on energy. 

As known, the slope parameter, B — d ln{da/dt)/dt\t^Q, is increasing function. It is mainly connected 
with irdq/ sinh{TTdq), and with the parameter d. So, the parameter d must be increasing function also. 
The same regularities can be seen for pp-interactions. 

The fitting results for R and d are presented in Fig. 5. They show that R for pp-interactions decreases 
with the energy growth starting from low energy, reaches a minimum at ^/s ^ 30 GeV, and continues the 
growth at higher energies. R for pp-interactions in the studied energy range is practically constant. 

The energy dependence of d is more complicated. For pp-interactions in the considered energy range 
it is the increasing functions, d for pp-scattering has an interesting irregularity at very low energies. At 
high energies it reaches a constant value. 
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Figure 5: The fitting results for the parameters R and d. The black points present the results for 
pp-interactions, the open points - for pp-interactions. 

For future applications we approximate the dependencies as: 

^ 0.07 + 0.05 lns+^ + i|:^, (/m), (14) 

= 0.07 + 0.05 In s+i^ + -|^, (/m), (15) 
0.5 

dpp^dpp^ 0.36 J^. ifm). (16) 

The dependencies are shown in Fig. 5 by the solid and dashed lines. The asymptotical part, 0.07 + 
0.051ns, is presented by the dotted line. 

The fitting results for the parameters A and oi are shown in Fig. 6. As seen, the parameter A 
fluctuates within ±30 % at low energies. It is close to unity at high energies, and there is a defined energy 
dependence of the parameter. Thus, we believe that the parameter does not refiect only uncertainty of 
the experimental data normalization, but it contains some information on physics of the processes. 

If the parameter A is below unity, it points out on a possible influence of the inelastic shadowing on 
the elastic scattering due to the processes of excitations and deexcitations of low mass difractive states 
during the scattering. A value of the parameter above the unity can be interpreted as a presence of 
additional processes like 7r-meson exchange, annihilation and so on which are not taken into account 
directly. The value above unity can violate the unitarity requirement according to which \j{b)\ must be 
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below unity. If yl < 1 + e^^l'^ there is no problem with the unitariry, but it will mean that the amplitude 
reaches the black disk limit in the central interactions. If A > 1 + the simplest solution can be an 

application of any unitarisation scheme. We are going to study the subject in the future. 

We show in Fig. 6a the function 1 + e^^l'^ for pp-interactions by the dashed line. It seems to us that 
taking into account the fluctuation of the fitting results for the pp-interactions at high energies we can 
assume that the black disk limit is reached in the pp-collisions in the central region. 
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Figure 6: The fitting results for the parameters A and a\. The black points present the results for 
pp-interactions, the open points - for pp-interactions. 

As seen in Fig. 6, energy dependence of the parameter a\ for pp-interaction is rather regular. We 
cannot say this for pp-interactions. 

For future applications we approximate the dependencies of the parameters A and a\ as: 

^pp = 1 + e-^/^ = 1 + e-^^l^ ~ 0.6/s°-25, (17) 

ai, pp = 0.5/s, ai, PP -0.1/s. (18) 

Using the expressions [T^] - HH] we obtain a good description of the pp- and pp-interaction data (see 
Appendixes). Especially, we have for pp-data y^jNoF = 4866/1489 ~ 3.26 at Piah > 9.9 GeV/c, and for 
pp-data ~ /NoF = 3620/675 ~ 5.36 at Piab > 8 GeV/c having only one fitting parameter - A. 



2 Description of the Totem data 

The expressions [T^ HM [T51 predict for the LHC energies the following values of the parameters: 



yi(GeV) 


R (fm) 


d (fm) 


P 


(Ttot (mb) 


aei (mb) 


900 


0.770 


0.359 


0.1316 


71.5 


16.9 


7000 


0.963 


0.360 


0.1346 


90.9 


22.7 


10000 


0.997 


0.360 


0.1347 


94.8 


24.0 


14000 


1.030 


0.360 


0.1348 


98.8 


25.3 



Table 1: Parameters and estimated value for LHC 



Using them we calculate total and elastic cross sections, atot and aei- 
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(19) 



CTe; ~TT {R^ - 19M/10 + 2n^(f/7). (20) 

The Totem collaboration [33] published the total and elastic cross sections which are atot — 98.3 ± 
Q 2stat_^2.8^2'' (mb) and a^i = 24.8±0.2''*°*±1.2^2^^ (mb). They are above our predictions. To understand 
the difference, we calculate differential elastic scattering cross section according to our approach, and the 
cross section using the simple gaussian parameterization and the value of the slope parameter given by 
the CoUaboration, B = 20.1 ± 0.2"*"* ± 0.3"^" (GeV^^). They are presented in Fig. 7 by the sohd and 
dashed lines, respectively. 




|t| (GeV') 

Figure 7: Differential cross section of pp elastic scat- 
tering at 7 TeV. The points are the experimental 
data [Ml US] • The lines are our calculation. 

As seen, first of all, the Collaboration fitted the differential cross section at \t\ < 0.35 (GeV~^) to 
obtain the slope and the cross sections (dashed line). Our prediction (sohd line) catches the points at \t\ > 
0.15 (GeV~^, especially, in the region of the minimum. At |t| below 0.15 (GeV~^) the prediction deviates 
regularly from the corresponding dataQ. The data are above our curve. Maybe additional expansion 
terms are needed to be included in Eq. [11] The can give corrections at small 

Our description of the data [24] [23] in the whole measured values of t is presented in Fig. 8. 




O.O 0.5 1.0 1.5 2.0 2.5 



|t| (GeV^) 

Figure 8: Differential cross section oipp elastic scattering at 7 TeV. The points are the experimental data 
[24l [23] . The line is our calculation. 

*We could not be able to digitize quite well the experimental points presented in 23 . Thus experimental errors are not 
shown. We cannot guarantee that exactness of the shown points is sufficiently high. 
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One can see that the forward scattering data are reproduced quite well. The dip is filled rather well. 
Its position is right. At the same time, the calculations deviate from the data starting from lower values 
of |i| than it was at other energies. The high of the second diffractional maximum is mainly determined 
by the parameter d. The slope of the forward part of the spectra is connected with the parameter also. 
At chosen value of the parameter we overestimate a little bit the high of the maximum. We can describe 
better the forward part of the date varying d in its accuracy limits making worse the description of the 
dip region, and vice- versa. We expect that an exactness of the parameter determination will be improved 
when the Totem collaboration will publish final data. 

We have to note that an accuracy of the parameters entering in Eqs. [2] -[111 HHis equal to ±5 %. 
We expect the same accuracy for the calculated differential cross section. The accuracy can be improved 
when new Totem data at other energies will be appeared. 

In order to understand a quality of the calculations, let us compare our calculations with predictions 
of other models [5S] presented by the Totem collaboration in the paper [M]. For this, we show the model 
predictions in the dip region and in the region of large |<| in Fig. 9. Because we could not be able to 
take experimental errors in [24_, we plotted the points without errors. Though, the systematic errors are 
rather large for a correct discrimination of the model, one can see that only our approach gives results 
that are quite closed to the data at small angles and in the dip region. The high of the second maximum 
is reproduced also in the approach. But instead of other models we predict too fast decreasing of the 
cross section above the diffraction maximum. The other models predict much slower decreasing of the 
cross sections. Here we have to note, that the models were tuned using much less set of experimental 
data. Additional to this, they included, directly or indirectly, the high momentum transfer scattering. 




|t| (GeV') |t| (GeV') 

Figure 9: Comparison of model predictions presented in [M] with the Totem experimental data. The 
solid (black) lines show our calculations. 



3 Description of the high momentum part of the data 

The nature of processes with high momentum transfer are debated until now. It is commonly accepted 
that they can be described in QCD. There are a lot of publications on the subject. Instead of analyzing 
all of them in order to select a reliable one we turn to experimental data. In Fig. 10 we present some 
experimental data on elastic scattering at large 

As seen, at \t\ > 2 (GeV/c)^ aU the cross sections have the same shape at y/s > 10 GeV. At the 
projectile momentum below 200 GeV/c they have strong energy dependence. To reproduce the high 
energy behavior of the cross sections we add to the imaginary part of the amplitude (|lip a "hard" 
scattering amplitude: 

ImUardiq) = -0.05 [1 + tanh{Rq - ^■^)]-(Y+W^YT)^- ^^^^ 

The sign "-" is needed to increase the cross sections in the second maximum where Ji{Rq) is negative. 
The hyperbolic tangent imitates a smooth transition from soft to hard scatterings. According to Fig. 10 
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the border of the hard processes slowly moves to small momentum transfer with an energy growth. It can 
be if the border is connected with the radius of the soft interactions. Thus, we assume that the tangent 
argument is Rq. The value "5.5" gives the exact position of the border. The last factor in Eq. [5T]is the 
proton form-factor in a tuned power. All values in (|21[) are sampled only in order to reproduce the cross 
section behavior qualitative at > 10 GeV. With all of these we have a description of the Totem data 
presented in Fig. 11. 




|t| (GeV') |t| (GeV') |t| (GeV') 



Figure 10: pp elastic scattering data at large momentum transfer. Points are experimental data 

mmmmmmm- 

Of course, our parameterization of the high momentum part is not perfect one. But at least, it 
describes the previous experimental data, and we cannot simple disregard it. The behavior of the pre- 
dictions in Fig. 11 is explained by the variation of R. As energy increases, R is increased also, and the 
yield of the soft part in the high |t| region decreases, the dip is shifted to the lower |t|, and the second 
maximum increases. If we are right, the future measurement of the Totem collaboration can show this. 
The measurement will give us more information about interplay of the soft and hard interactions. 




0.0 0.5 1.0 1.5 2.0 2.5 

|t| (GeV^) 

Figure 11: Description of the Totem data and our predictions for other LHC energies. 
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4 Concluding remarks 

4.1 Total and elastic cross sections 

According to Exps. ^ and © the total cross section is given by the expression: 

atot ^ 27r A{R^ + TT^(f/S). (22) 
For the elastic cross section we have the following expression an exactness of which is about few percent. 

<7ei c^n {R^ -19 Rd/ 20 + n^(f/ 7). (23) 
They together with Exps. (ITil) - ([T7)) provide one with a good parameterization of the cross sections. 

4.2 Eikonal representation 

The profile-function (jlOp can be represented as an eikonal one: 

1 _ e" H(e'"''+e''''')/{e''^''-l)] _ (^24) 



lib) = 1 



I + g(h-i?,)/d 



Here we use the assumption that A = 1 + e^^/''-. 

It can be used in the Quark-Gluon String Model 'S^ for a calculations of string multiplicity distribu- 
tions. 

P„ = i f^ {^H{e'''" + e^")/{e'"-lW ^-2I„[(e«/^+eV^)/(eV^-l)] ,25) 

C J n! 

4.3 Application of USESD in calculations of hadron-nucleus and nucleus- 
nucleus properties. 

An amplitude for an elastic scattering of an nucleus containing B baryons on a target nucleus with mass 
number A is given as [28] : 

FBAiq)^^ J d^e^^^ i^l-f[J[[l-^{Un-s,)]^\^B\'\^A\' {flf*^ (P/''') ^^^^ 

bPBA{b) Mqb)db, 

where 7 is an amplitude of an elastic nucleon-nucleon scattering in the impact parameter representation, 
averaged over the spin and isospin degrees of freedom, 

lib) -^Jd'q e^'^fiq). 

^A i^B) is the wave function of the target (projectile) nucleus in the ground state. Taking the origin of 
the coordinate system to coincide with the center of the nucleus, the nucleon coordinates ({r^}, {is}) 
are decomposed into longitudinal ({^i}) and transverse ({sj}, {fi}) components. The z-axis is directed 
along the projectile momentum, b is the impact parameter vector orthogonal to the momentum. Pba ib) 
is the profile function and Jq is the Bessel function of zero order. 
Quite often 7 is parameterized a^: 

7(&) = ^^''^^ e-''/'f>, (27) 

where cr^^ is the total cross section of the nucleon-nucleon interactions, p is the ratio of the real to 
imaginary parts of the NN elastic scattering amplitude at zero momentum transfer, and /? is a slope 
parameter of the NN differential elastic scattering cross section. Then 

fiq)-{^<Nii-^P)e-^'"^^ (28) 



^The amplitude must be corrected at low energies in order to take into account the unitarity requirement (Re7(0) < 1) 
and a restriction of the phase space. 
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j3 can be found as: 

P - + \P?)im ^ ^NN 0.3897). 

Here, cr^jy elastic cross section and 0.3897 is a coefficient required in order to express /3 in 

units of (GeT^/c)"^, if the cross sections are given in millibarns. 
The squared modulus of the wave function is usually written as: 

A A 

\^A\'^S{J2r./A) l[pA{r.). (29) 

i=l i=l 

PA coincides with the one-particle density of the nucleus if one neglects the center-of-mass correlation 
connected with the S function. 

We have used Exp. ([27)) at the calculation of the differential elastic scattering cross section at ^/s — 7 
TeV, results of which are shown in Fig. 7 as the dashed line (Gaussian param.). It is obvious, that it can 
describe the cross sections only at small value of It would be better to use Exp. (flUl) for more exact 
calculations as nucleon-nucleon interaction properties, as well as baryon-nucleus and nucleus-nucleus ones. 

The task is very actual for calculations of wounded nucleon multiplicities in nucleus-nucleus interac- 
tions at high energies. 



4.4 Comparison of USESD with other approaches 

Many years ago T.-Y. Cheng, S.-Y. Chu and A.W. Hendry successfully used 2-dimensional Fermi- 
function to describe pp elastic scattering at all angles from 3 to 24 GeV/c. They analyzed polarized 
proton scattering by proton. Most probably is that they used numerical integration. We have used 
analytical expressions, and we have considered only unpolarized proton scattering. 

In 2000 M. Kawasaki, T. Maehara and M. Yonezawa analyzing the general structure of the elastic 
scattering amplitude in the impact parameter representation proposed the following expression for the 
amplitude in the momentum representation [30) : 



/m f{q) cx 



Rq pq 

where is a parameter. They proposed also a concrete form of the dumping functions: 
<i>o{R,q) = i?e r [l + i-] , '^i{R,q) = -Im T ( 1 + i- 



(30) 



They were trying to fit the high energy pp- and pp-data [SU |32] but results were not impressive. 
Especially, in the paper [32] the authors fitted the differential cross sections of pp elastic scattering at 
y/s > 23.5 GeV and the pp-data at ^/s > 546 GeV. Only small momentum transfer region, 0.02 < \t\ < 0.2 
(GeV/c)^, was included in the fit. 

In our paper we have considered much larger energy region and a larger region of the momentum 
transfer. We assume our results as promising ones. Thus a choice of the dumping functions is very 
important for a correct reproduction of experimental data. 

P. Gauron, B. Nicolescu and E. Leader proposed in the paper [3^ the following expressions for 
asymptotic parts of (anti) proton-proton elastic scattering amplitude: 

-F+{s,t) = Fi ln\s)^^^4r^ e''^* +F2 ln{s)MR+f) e''^* + (31) 

iS R+T 

Fs [MR+f) - R+f JiiR+f)] e''^*, 
-F_{s,t) = Oi ln^{s)^^^4r^ e^^* + O2 lnis)cosiR-f) e^^* +O3 6^3', (32) 

S R+T 

where s = ^e~*'^"; r = (— j ln{s); sq = to = 1 GeV^; F^, Oi, hi, R+ and i?_ are constants. 

Fpp — F+ -f F_ , Fpp — F+ — F-. 
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The appearance of the exponents in the expressions is connected with the simplest assumption about 
the functional form of the residue functions. If one replaces the exponents by nd^/^ / sinh{TTd^/^) then 
the leading term of will be coincided with the first term of our expression (jlip . Thus our approach 
corresponds to the approach of the authors at defined assumptions on the residue functions. But instead 
of the authors we did not assume any energy dependence of our parameter R. We obtained it at the 
fitting of the experimental data. 

If we suppose that R = Rq + AR where AR cx s~°''^-, and expand the first term of Eq. [11] we will 
have yields of "non-dominant Regge pole contributions" into the scattering amplitude. 

Very often in a Regge-like analysisH the elastic scattering amplitude of a process 1 + 2— >• 1 + 2 is 
represented as: 

lib) = T^{l-e--^'^), (33) 

where 

^ Rl + R'^ + a'^' ' > 

rj = l+icotan(TTap{0)/2) is a signature factor; Ci (C2) - shower enhancement coefficient in the interaction 
vertex of the first (second) particle with a reggeon/pomeron; ap{0) = 1 + A is an intercept of the 
reggeon/pomeron; ^ = In S12 is a logarithm of CMS energy squared; C = £, ~ *7r/2. It is assumed that 
Regge trajectories are linear, a{t) = 1 + A + ^'t. Non-linear trajectories were considered in [33]. It is 
assumed also that the residue functions have the gaussian shape - gi{t) = gi{0)exp{R^ t), for simplicity. 

It is complicated to find a correspondence between the eikonal [34l and our eikonal [24l Though, the 
structure of the eikonal [34] is rather simple - it is a product of a function of 6, and a factor strongly 
dependented on the energy - e^^. For our eikonal at large b and Ci = C2 = 1, we have: 

Lu{b) = ln[(e^/'^ + e''/'')/(e''/'^ - l)]|b^oo ^ e-''/'*[e^/'^ + 1]. 

Thus the energy dependence of the eikonal is determined by e-'^'^'^ ~ gO.i/d ^ 5O.28 jjgj-g take into 
account that R cx 0.11ns, and d cx 0.36. So, an effective intercept in our model is 1.28. It is in a 
correspondence with results of papers |35l 136] . 

Summing up, we can say that our model is in the main stream of phenomenological analysis of the 
elastic scattering data. The model assumes the defined choice of the dumping functions, or the residue 
functions. A correct form of the functions is very important for high energy phenomenology. 

The authors are thankful to the Geant4 hadronic working group for interest in the work. 
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Appendix A: Comparison of experimental 
ization 
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Figure 12: The points are the experimental data 
by R.M. Edelstein et al., Phys. Rev. D5 (1972) 
1073. 
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Figure 14: The points are the experimental data 
by G.W. Brandenburg et al., Phys. Lett. 58B 
(1975) 367. 



a on pp-interactions with USESD parameter- 
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Figure 13: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 
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Figure 15: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425. 
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Figure 16: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 15 (1965) 
45. 
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Figure 17: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 
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Figure 18: The points are the experimental data 
by D. Harting, Nuovo Cimento 38 (1965) 60. 
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Figure 19: The points are the experimental data 
by K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425. 
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Figure 20: The points are the experimental data Figure 21: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) by K.J. Foley et al, Phys. Rev. Lett. 11 (1963) 
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Figure 22: The points are the experimental data 
by R.M. Edelstein et al, Phys. Rev. D5 (1972) 
1073. 
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Figure 24: The points are the experimental data 
by D. Harting, Nuovo Cimento 38 (1965) 60. 
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Figure 26: The points are the experimental data 
by .K.J. Foley et al., Phys. Rev. Lett. 11 (1963) 
425 
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Figure 23: The points are the experimental data 
by .K.J. Foley et al, Phys. Rev. Lett. 11 (1963) 
425 
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Figure 25: The points are the experimental data 
by J.V. Allaby et al., Phys. Lett. B28 (1968) 67. 
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Figure 27: The points are the experimental data 
by R.M. Edelstein et al., Phys. Rev. D5 (1972) 
1073. 
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Figure 28: The points are the experimental data 
by J.V. Allaby et al., Phys. Lett. B28 (1968) 67. 
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Figure 30: The points are the experimental data 
by R.M. Edelstein et al, Phys. Rev. D5 (1972) 
1073. 
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Figure 32: The points are the experimental data 
by C. Bruneton et al., Nucl. Phys. B124 (1977) 
391. 
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Figure 29: The points are the experimental data 
by J.V. Allaby et al., Nucl. Phys. B52 (1973) 
316. 



> 

0) 

O 
E 

"B 
•a 



10' 



10" 



10' 



10-' 



10-^ 



lo-' 



P, =35 GeV/c 

lab 



PP 



0.0 



0.5 



1.0 



1.5 



2.0 



|t| (GeV^) 



Figure 31: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 33: The points are the experimental data 
by Z. Asad et al., Nucl. Phys. B255 (1984) 273; 

C. W. Akerlof et al., Phys. Rev. D14 (1976) 2864; 

D. S. Ayres et al., Phys. Rev. D15 (1977) 3105. 
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Figure 34: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 35: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 36: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 
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Figure 37: The points are the experimental data 
by C.W. Akerlof et al., Phys. Rev. D14 (1976) 
2864; R. Rubinstein et al., Phys. Rev.DSO (1984) 
1413. 
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Figure 38: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105; R. Rusack et al., Phys. Rev. Lett. 41 
(1978) 1632. 
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Figure 39: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
191632. 
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Figure 40: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 
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Figure 41: The points are the experimental data 
by D. Brick et al., Phys. Rev. D25 (1982) 2794. 
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Figure 42: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 



Figure 43: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 44: The points are the experimental data 
by D.S. Ayres et al., Phys. Rev. D15 (1977) 
3105. 



Figure 45: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 46: The points are the experimental data 
by C.W. Akerlof et al., Phys. Rev. D14 (1976) 
2864; R. Rubinstein et al., Phys. Rev.DSO (1984) 
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Figure 48: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 50: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 47: The points are the experimental data 
by A. Schiz et al., Phys. Rev. D24 (1981) 26; G. 
Fidecaro et al., Nucl. Phys. B173 (1980) 513. 
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Figure 49: The points are the experimental data 
by R. Rusack et al., Phys. Rev. Lett. 41 (1978) 
1632. 
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Figure 51: The points are the experimental data 
by M.G. Albrow et al., Nucl. Phys. B108 (1976) 
1. 
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Figure 52: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 53: The points are the experimental data 
by M.G. Albrow et al., Nucl. Phys. B108 (1976) 
1. 
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Figure 54: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 55: The points are the experimental data 
by A. Breakstone et al. Nucl. Phys. B248 (1984) 
253. 
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Figure 56: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 



Figure 57: The points are the experimental data 
by E. Nagy et al., Nucl. Phys. B150 (1979) 221. 
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Figure 58: The points are the experimental data 
by A. Breakstone et al. Nucl. Phys. B248 (1984) 
253; Phys. Rev. Lett. 54 (1985) 2180. 
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Figure 59: The points are the experimental data 
by S. Erhan et al., Phys. Lett. B152 (1985) 131; 
J.C.M. Armitageetal.,Nucl. Phys. B132 (1978) 
365. 
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Figure 62: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 



Figure 63: The points are the experimental data 
by U. Amaldi and K.R. Schubert, Nucl. Phys. 
B166 (1980) 301. 
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Figure 64: The points are the experimental data 
by J.S. Russ et al., Phys. Rev. D15 (1977) 3139. 
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Figure 66: The points are the experimental data 
by A. Berglund et al., Nucl. Phys. B176 1980) 
346. 
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Figure 67: The points are the experimental data 
by G. Brandenburg et al., Phys. Lett. 58B 
(1975) 367. 
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Figure 68: The points are the experimental data 
by K.J. D. Birnbaum et al., Phys. Rev. Lett. 23 
(1969) 663. 
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Figure 70: The points are the experimental data 
by Yu.M. Antipov et al., Nucl. Phys. B57 (1973) 
333. 
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Figure 72: The points are the experimental data 
by Yu.M. Antipov et al, Nucl. Phys. B57 (1973) 
333. 
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Figure 69: The points are the experimental data 
by B. Batyunya et al. Yad. Fiz. 44 (1986) 1489. 
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Figure 71: The points are the experimental data 
by M.Y. Bogolyubsky et al, Yad. Fiz. 41 (1985) 
1210. 
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Figure 73: The points are the experimental data 
by R.L. Cool et al., Phys. Rev. D24 (1981) 2821. 
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Figure 74: The points are the experimental data 
by A. Breakstone et al., Nucl. Phys. B248 (1984) 
253, Phys. Rev. Lett. 54 (1985) 2180. 
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Figure 76: The points are the experimental data 
by A. Breakstone et al., Nucl. Phys. B248 (1984) 
253, Phys. Rev. Lett. 54 (1985) 2180. 
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Figure 78: The points are the experimental data 
by M. Bozzo et al, Phys. Lett. B155 (1985) 197, 
Phys. Lett. B147 (1984) 385. 
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Figure 75: The points are the experimental data 
by A. Breakstone et al., Nucl. Phys. B248 (1984) 
253, Phys. Rev. Lett. 54 (1985) 2180. 
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Figure 77: The points are the experimental data 
by G. Arnison et al., Phys. Lett. B128 (1983) 
336. 
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Figure 79: The points are the experimental data 
by D. Bernard et al, Phys. Lett. B198 (1987) 
583, Phys. Lett. B171 (1986) 142. 
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Figure 80: The points are the experimental data 
by N. Amos et al., Nucl. Phys. B262 (1985) 689, 
Phys. Lett. B247 (1990) 127. 



Figure 81: The points are the experimental data 
by the DO Collaboration - DO Note 6056-CONF. 
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Appendix C: Quality of the data fitting 

presented below is obtained with values of the parameters given by Eqs. [12] ^ HH [181 ^ is nmnber 
experimental points. Np - the number of points included in the fitting. 



Table 2: Fitting results for pp- interactions 



Plab 

(GeV/c) 


1 — 

(GeV) 


N 


Np 


/NoF 


A 


R 
(fm) 


d 
(fm) 


P 


(fm^) 


9.9 


4.523 


28 


28 


0.560 


1.180± 0.018 


0.547 





249 


-0.323 


0.0244 


10 


4.543 


24 


21 


8.66 





746 


± 





016 


0.547 





250 


-0.322 


0.0242 


10.4 


4.625 


97 


97 


2.56 


1 


244 


± 





002 


0.545 





252 


-0.319 


0.0234 


10.8 


4.705 


13 


13 


2.70 


1 


177 


± 





012 


0.544 





254 


-0.315 


0.0226 


10.9 


4.733 


9 


9 


0.72 


1 


108 


± 





014 


0.544 





254 


-0.313 


0.0223 


12 


4.939 


25 


18 


8.48 





845 


± 





019 


0.540 





259 


-0.303 


0.0205 


12.4 


5.014 


22 


21 


0.72 


1 


065 


± 





012 


0.539 





260 


-0.299 


0.0199 


12.8 


5.088 


13 


13 


2.96 


1 


180 


± 





013 


0.538 





262 


-0.296 


0.0193 


14.2 


5.34 


22 


12 


4.32 





818 


± 





021 


0.535 





266 


-0.283 


0.0175 


14.8 


5.444 


12 


12 


1.98 


1 


159 


± 





014 


0.533 





268 


-0.278 


0.0169 


15.1 


5.496 


32 


25 


0.22 


1 


200 


± 





019 


0.533 





269 


-0.276 


0.0166 


16.7 


5.762 


12 


12 


1.22 


1 


154 


± 





014 


0.530 





273 


-0.263 


0.0151 


18.4 


6.033 


20 


18 


1.32 


1 


076 


± 





016 


0.527 





277 


-0.251 


0.0137 


19.2 


6.156 


27 


19 


14.4 





899 


± 





Oil 


0.526 





279 


-0.245 


0.0132 


19.6 


6.217 


11 


11 


1.82 


1 


132 


± 





014 


0.526 





280 


-0.243 


0.0129 


20 


6.277 


30 


28 


2.92 


1 


043 


± 





016 


0.525 





280 


-0.240 


0.0127 


21.12 


6.442 


17 


15 


12.04 





946 


± 





Oil 


0.524 





282 


-0.233 


0.0121 


24 


6.849 


35 


20 


5.62 


1 


000 


± 





015 


0.521 





287 


-0.217 


0.0107 


29.7 


7.591 


29 


25 


1.60 


1 


060 


± 





018 


0.517 





294 


-0.190 


0.0087 


35 


8.221 


21 


19 


1.30 


1 


120 


± 





022 


0.514 





299 


-0.171 


0.0074 


44.5 


9.243 


27 


25 


7.28 


1 


093 


± 





004 


0.511 





306 


-0.143 


0.0059 


50 


9.787 


33 


27 


11.06 


1 


106 


± 





006 


0.510 





309 


-0.130 


0.0052 


55.4 


10.29 


20 


17 


1.38 


1 


047 


± 





022 


0.510 





311 


-0.119 


0.0047 


65 


11.13 


20 


17 


1.14 


1 


134 


± 





024 


0.509 





315 


-0.102 


0.0040 


70 


11.55 


17 


17 


2.26 


1 


136 


± 





006 


0.509 





317 


-0.095 


0.0037 


100 


13.78 


81 


69 


1.22 


1 


147 


± 





007 


0.509 





324 


-0.062 


0.0026 


100 


13.78 


37 


34 


1.62 


1 


125 


± 





006 


0.509 





324 


-0.062 


0.0026 


130 


15.69 


11 


11 


0.34 


1 


137 


± 





030 


0.510 





328 


-0.040 


0.0020 


140 


16.28 


19 


19 


3.56 


1 


124 


± 





006 


0.510 





329 


-0.034 


0.0019 


150 


16.85 


31 


26 


0.60 


1 


107 


± 





021 


0.511 





330 


-0.029 


0.0018 


170 


17.93 


29 


23 


0.60 


1 


154 


± 





024 


0.512 





332 


-0.020 


0.0016 


175 


18.19 


15 


15 


2.82 


1 


139 


± 





007 


0.512 





333 


-0.018 


0.0015 


190 


18.95 


30 


24 


1.18 


1 


098 


± 





023 


0.513 





334 


-0.012 


0.0014 


200 


19.44 


90 


58 


1.44 


1 


104 


± 





007 


0.513 





334 


-0.009 


0.0013 


200 


19.44 


167 


156 


2.48 


1 


080 


± 





002 


0.513 





334 


-0.009 


0.0013 


210 


19.91 


29 


23 


0.72 


1 
i 


uyo 


1 

± 


u 


nop; 


0.514 





335 


-0.006 


0.0013 


230 


20.84 


19 


16 


2.34 


1 


064 


± 





030 


0.515 





336 


0.000 


0.0012 


250 


21.72 


17 


14 


0.68 


1 


167 


± 





033 


0.516 





337 


0.005 


0.0011 










n OP. 
U./D 


1 


186 


± 





024 


U.oio 





339 


U.Ui4 


u.uuuy 


293.5 


23.53 


133 


97 


5.56 


1 


067 


± 





003 


0.518 





339 


0.015 


0.0009 


497 


30.6 


15 


15 


0.56 


1 


249 


± 





027 


0.527 





344 


0.041 


0.0005 


501.3 


30.73 


124 


88 


2.32 


1 


135 


± 





003 


0.527 





344 


0.042 


0.0005 


512 


31.05 


24 


24 


0.48 


1 


100 


± 





015 


0.528 





344 


0.043 


0.0005 


1064 


44.74 


65 


26 


3.80 


1 


032 


± 





008 


0.544 





349 


0.070 


0.0002 


1486 


52.87 


63 


19 


3.48 





982 


± 





014 


0.553 





351 


0.080 


0.0002 


1497 


53.07 


55 


41 


1.20 


1 


027 


± 





013 


0.553 





351 


0.080 


0.0002 


1497 


53.07 


27 


24 


3.82 


1 


143 


± 





009 


0.553 





351 


0.080 


0.0002 


2048 


62.06 


23 


23 


0.58 


1 


094 


± 





013 


0.562 





352 


0.087 


0.0001 


2077 


62.5 


74 


37 


4.20 


1 


000 


± 





014 


0.563 





352 


0.088 


0.0001 


2081 


62.56 


49 


49 


13.6 





996 


± 





001 


0.563 





352 


0.088 


0.0001 
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Table 3: Fitting results for pp-interactions 



Plab 

(GeV/c) 


(GeV) 


N 


Nf 


x'/NoF 


A 


R 
(fm) 


d 
(fm) 


P 


ai 
(fm^) 


8.00 


4.111 


41 


41 


4.42 





833 


± 





0008 


0.895 


0.238 


-0.320 


-0 


005916 


9.71 


4.483 


10 


5 


1.90 





618 


± 





0427 


0.855 


0.249 


-0.320 


-0 


004975 


10.1 


4.564 


35 


31 


5.04 





880 


± 





0042 


0.847 


0.250 


-0.320 


-0 


004801 


10.4 


4.625 


61 


61 


3.68 





881 


± 





0008 


0.842 


0.252 


-0.320 


-0 


004675 


16.0 


5.647 


23 


23 


2.38 





834 


± 





0023 


0.766 


0.272 


-0.265 


-0 


003136 


22.4 


6.626 


44 


44 


1.28 





910 


± 





0047 


0.716 


0.285 


-0.206 


-0 


002277 


25.2 


7.012 


33 


33 


2.00 





885 


± 





0040 


0.701 


0.289 


-0.187 


-0 


002034 


32.1 


7.882 


52 


52 


4.98 





893 


± 





0035 


0.673 


0.297 


-0.152 


-0 


001609 


40.1 


8.785 


33 


33 


1.08 





888 


± 





0068 


0.650 


0.303 


-0.122 


-0 


001296 


191.1 


19.1 


27 


23 


2.28 





974 


± 





0182 


0.562 


0.334 


0.016 


-2 


770E-4 


510.2 


31.0 


22 


22 


1.44 


1 


076 


± 





0080 


0.550 


0.344 


0.062 


-1 


041E-4 




53.0 


51 


44 


3.94 


1 


015 


± 





0053 


0.560 


0.351 


0.092 


-3 


560E-5 




62.3 


23 


23 


1.12 


1 


068 


± 





0063 


0.566 


0.352 


0.099 


-2 


576E-5 




540 


36 


36 


0.84 


1 


090 


± 





0062 


0.720 


0.359 


0.131 


-3 


429E-7 




546 


121 


121 


17.24 


1 


143 


± 





0018 


0.721 


0.359 


0.131 


-3 


354E-7 




630 


19 


15 


7.66 


1 


124 


± 





0147 


0.734 


0.359 


0.131 


-2 


520E-7 




1800 


51 


51 


1.04 


1 


085 


± 





0028 


0.830 


0.360 


0.134 


-3 


086E-8 




1960 


17 


17 


0.36 


1 


095 


± 





0323 


0.838 


0.360 


0.134 


-2 


603E-8 
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Appendix D: Attempt to fit the COSY data on pp- interact ions D. Albers et al.,Phys. Rev. 
Lett. 78 (1997) 1652. 

We represent the elastic scattering amplitude in the energy range, 1-3 GeV/c as 

Fpp^fi9)+,f{7r-9), 

where f{9) is given by Eq. [TTJ Below the lines are results of a fit with 4 free parameters - A, R, d and 
ai. 
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Figure 82: The points are the experimental data. Figure 83: The points are the experimental data. 
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Figure 84: The points are the experimental data. Figure 85: The points are the experimental data. 
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Figure 86: The points are the experimental data. 
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Figure 87: The points are the experimental data. 
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Figure 88: The points are the experimental data. 
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Figure 89: The points are the experimental data. 
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Figure 90: The points are the experimental data. Figure 91: The points are the experimental data. 
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